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f— *^ , The medium modification of vector meson properties and its effects on thermal 

f*^ i photon and dilepton emission rates from hot/dense hadronic matter are discussed. 

' The role of spectral changes of hadrons in explaining the WA98 photon data and 

the CERES/NA45 dilepton data with similar initial conditions is demonstrated. 
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1 Introduction 



A fundamental property of QCD with two light quark flavours is the approximate 
chiral symmetry which plays a very important role in the understanding of low energy 
' hadron physics. Since this symmetry is not manifest in the observed particle degen- 

^sO . eracies it is believed that chiral symmetry is spontaneously broken by the vacuum. 

Under conditions of high temperature or density one expects that the QCD vacuum 
will undergo substantial modifications leading to the restoration of chiral symmetry. 
The chiral condensate which has a non-trivial value in the broken phase acts as the 
order parameter of this transition. Though hadron masses can be expressed in terms 
of various condensates through the QCD sum rules for example, there are uncertain- 
ties connected with the approximations used in their evaluation. The change in the 
condensate in hot/dense environment is expected to show up in the spectral function 
of the vector mesons. Concomitant to these phenomena, quarks and gluons, which 
are normally bound within hadrons are expected to undergo a deconfinement tran- 
sition at sufficiently high temperature or baryon density leading to the formation of 
quark gluon plasma (QGP)p. 

Relativistic collisions of heavy ions provide a rich arena for the study of strongly 
interacting matter under extreme conditions in the laboratory. Data from the Su- 
per Proton Synchrotron (SPS) have already been analysed and the first results from 
the Relativistic Heavy Ion Collider (RHIC) have started to appear. Electromagnetic 
probes, viz., photons .and dileptons have long been recognized as the most direct 
probes of the collisionB. Owing to the nature of their interaction they undergo min- 
imal scatterings and are by far the best markers of the entire space-time evolution 
of the collision. Both photons and dileptons couple to hadrons through spin one 
(J p = 1~) mesons. The dilepton spectra in particular, exhibits a resonant struc- 
ture which, in the low mass regime includes the p and the uj mesons. It has been 
emphasized that the properties of the vector mesons undergo nontrivial medium mod- 
ifications in a hot and/or dense medium such as likely to be produced in relativistic 
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nuclear collisions. Consequently the spectral modifications of these vector mesons 
would be clearly revealed in the invariant mass spectra of the dilcptons through the 
shift of the resonant peaks. 

In the following we will discuss the role of medium modification of hadronic 
properties in explaining the photon and dilepton spectra from heavy ion collisions 
at the CERN SPS. Let us first consider the data obtained from Pb-Pb collisions 
at 158 A GeV reported by the WA98 Collaboration El which shows the transverse 
momentum spectra of excess photons above those coming from the decay of ir° and 
r\ mesons. These excess photons are believed to arise due to (a) the hard collisions 
of partons in the colliding nuclei and (b) a thermal medium - QGP or hot/dense 
hadronic matter. We have seenta that the above spectra can be explained using a 
reasonable set of parameters if the medium effects of the vector meson properties are 
consistently incorporated in the invariant rates as well as in the equation of state. 
It is worth recalling that the single photon data for S+Au collisions could also be 
explained in a similar scenario □. However, other attempts tQjjVxrjlain the photon 
spectra without invoking medium effects exist in the literaturet] Blj U liil 

Our aim now is to check whether the same scenario can explain the dilepton data 
obtained in Pb+Au collisions at SPS by the CERES/NA45 Collaboration 0. The 
possible sources in the low mass region are the dileptons coming from the decays 
of hadrons at freeze-out and from the in-medium propagation and decay of vector 
mesons. The data shows_a significant enhancement in the mass region 0.3 to 0.6 GeV 
which can be explainedE3 by a substantial negative shift of the mass of vector mesons 
(the p meson in particular) in the thermal medium. A large broadening of the p meson 
spectral function due to scattering off baryons can also explain this enhancementH3. 

In the following we will see that both the WA98 photon spectra as well as the 
CERES/NA45 dilepton spectra for higher multiplicities can be explained by either 
of the two scenarios of relativistic nuclear collision liim: (a) A +A ^QGP^Mixed 
Phase^Hadronic Phase or (b) A + A ^Hadronic Matter; with downward shift of 
vector meson masses and initial temperature ~ 200pMeV. The photon yield is seen to 
be insensitive to the broadening of vector mesons although the CERES dilepton 
data admits such a scenario as mentioned above. The effects of the thermal shift 
of the hadronic spectral functions on both photon and dilepton emission have been 
considered in Ref. 113 for an exhaustive set of models. An appreciable change in 
the space-time integrated yield of electromagnetic probes was observed for universal 
scaling and Quantum Hadrodynamic (QHD) model and we will consider only these 
in our discussion. 

This report is organized as follows. In the following section we will briefly outline 
the phenomenology of medium effects on the vector mesons in the thermal environ- 
ment which we have considered. In section 3 we will consider the static photon and 
dilepton rates due to different processes. Then in section 4 we will describe the dy- 
namics of space-time evolution followed by the results of our calculation in section 5. 
We will conclude with a summary and discussions in section 6. 
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2 Medium Effects 



A substantial amount of literature has been devoted to the temperature 

and/or density dependence of hadrons within various modcistataotata in this work 

we will consider medium modifications of vector mesons in two different scenarios : 
the universal scaling hypothesis and the Quantum Hadrodynamic (QHD) model. 

In the scaling hypothesis, the parametrizatioryaf in-medium quantities (denoted 
by *) at finite temperature, T and baryon densityt3, Ub is 

m v _ fv _ w o f, „ „ 71 b\ f, T 2 



, - =1-0.2^ 1-^ , (1) 

where V stands for vector mesons, fv is the coupling between the electromagnetic 
current and the vector meson field, ojq is the continuum threshold, T c is the critical 
temperature and n B is the baryon density of normal nuclear matter. Mass of the 
nucleon also varies with temperature according to Eq. ([j]) (pseudo scalar masses 
remain unchanged). i-The particular exponent A = 1/2 (1/6) is known as Nambu 
(Brown- Rho) scalinglI3. We will use A = 1/2 in our calculations. Note that there is 
no definite reason to believe that all the in-medium dynamical quantities are dictated 
by a single exponent A. This is the simplest possible ansatz (seeta for a discussion). 
The effective mass of a\ is estimated by usipa-Weinberg's sum rulescJ. 

In the Quantum Hadrodynamic model r 1 tr 2 l of nuclear matter the vector meson 
properties are modified due to coupling with nucleonic excitations. The nucleons 
interact through the exchange of scalar a and the vector w mesons and their mass is 
modified due to the scalar condensate. This is evaluated in the Relativistic Hartree 
Approximation (RHA). Coupling with these modified nucleonic excitations induce 
changes in the p and ui meson masses. This modification is contained in the meson self 
energy which appears in the Dyson- Schwinger equation for the effective propagator 
in the medium. The interaction vertices are provided by the Lagrangian 



-VNN 



9vnn (n^NY* - -^L-Na^r a Nd»V a ^ , (2) 



where — {w^jp 1 "}, Mn is the free nucleon mass, N is the nucleon field and 
r a = {1, r}, f being the Pauli matrices. The real part of the veplatjaeson self-energy 
due to NN polarization is responsible for the mass shift (seeHt3n3'E3). The p — tt 
interaction is given by 



= " Wm ' (* x + «9%„ (t? x flj • (tt x (F). (3) 



This makes the dominant contribution to the imaginary part of the p self-energy 
though the pole shift (real part) is negligibly small. We have evaluated the spectral 
function of p due to its coupling with the pion and nucleons (Eqs. (|^) and (^)) in the 
thermal bath for finite three momentum at non-zero temperature and density. 
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3 Thermal Photon and Dilepton Emission Rates 



In this section, we briefly recapitulate the main equations relevant for evaluating pho- 
ton and dilepton emission from a thermal source. The emission rate of real photons 
with four momentum p M = (E,p) can be expressed in terms .of the imaginary part of 
the retarded electromagnetic current correlation function ascl 

E Sp- = ~ J^y r ImW * M fBE{E ' T) ' (4) 

where 

W* (p) = ij d 4 x e** 6(x ) ([J e ™(x), JT(P)]h > (5) 

J^ m being the electromagnetic current of quarks or hadrons. All the information 
about the thermal medium which produces the photons or dileptons resides in the 
current correlator W^ v . The lowest order processes contributing to thermal pho- 
ton emission from quark gluon plasma .are the QCD Compton and annihilation pro- 



cesses Ell. It has been shown recently E3 that the two-loop contribution leading to 
bremsstrahlung and qq annihilation with scattering is of the same order as the lowest 
order processes. The total rate of emission per unit four-volume at temperature T is 
given by 



dR 5 aa s 



In 



2.912 E 



(6) 



where Jt ^ 4.45 and Jl — -4.26. 

Considering the hadronic matter to be made up of n, p, co, rj and a\ mesons, the 
imaginary part of W^ v gives the amplitudes of a number of hadronic processes. We 
have considered the reactions n p — ► 717, tttt — > pj, tttt — > rjj, wrj — > and the 
decays p — > n n 7 and u> — > n 7 with dressed vector propagators to estimate photon 
emission from hadronic matter. 

To obtain the dilepton emission rate, the additional part to be included in Eq. (Q) 
corresponds to the free propagation of the virtual photon and its subsequent decay 
to lepton pairs. This gives, 

^=-T^7^ W >)( 1 + f ^ L ) \ll-^fBB[E,T) (7) 

where m; is the lepton mass. 

In the QGP, the dominant channel for dilepton productions is quark-antiquark 
annihilation. The rate for this process is obtained from the lowest order diagram 
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contributing to the current correlator W^ v and is obtained as, 



dR \ 2 a 2 




BE 



(E,T) 



(8) 



Now let us consider dilepton production in the hadronic medium. In order to 
obtain the rate of dilepton emission from hadronic matter (p/tu — > from Eq. (fij) 

the electromagnetic current correlator is expressed in terms of the effective propagator 
of the vector particle in the thermal medium using vector meson dominance (VMD) 
so as to obtain 
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(9) 



where Imlly is the transverse(longitudinal) part of the retarded self-energy of 
vector mesons arising out of interaction with excitations in the medium. The dilepton 
rate from pion annihilation for example, can be obtained from p self-energy due to 



7r polarisation in the medium. In the approximation IH = II 



= Hp we have 
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(10) 



where m* 2 = m 2 — Rellp and p 2 = M 2 , the invariant mass of the lepton pair. This 
is proportional to the familiar electromagnetic form factor of the pion. 



4 Space-Time Evolution 

The observed spectrum originating from an expanding hadronic matter is obtained 
by convoluting the static rates given above with the expansion dynamics. This is 
done using (3+1) dimensional relativistic hydrodynamics assuming boost invariance 
in the longitudinal direction and cylindrical symmetry in the transverse plane. The 
initial energy density and transverse velocity profiles are taken as 

£ fo' r > %(r-«S/« + i ; Vr=V °{~k) 

with a = 1. The other essential input one requires is the equation of state (EOS) 
which provides the cooling law. For the QGP sector we use the bag model equation of 
state with two flavour degrees of freedom. The hadronic phase is taken to consist of tt, 
p, u>, r\ and a\ mesons and nucleons. The medium effects enter through the effective 
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Figure 1: Total photon yield in Pb + Pb collisions at 158 A GeV at CERN-SPS. The system is 
formed in the QGP phase with initial temperature T{ = 190 MeV. 



masses in the expressions for energy density and pressure. The entropy density is 
then parametrized as, 

s H = = 4a cff (T) T 3 = 4^g cS (m*, T)T 3 (12) 

where g g is the effective statistical degeneracy. Thus, we can visualize the finite 
mass of the hadrons having an effective degeneracy g e g(m* ,T). The velocity of 
sound which plays a significant role in the expansion, also becomes a function of 
T and differs substantially from its value corresponding to an ideal gas (1/V3). The 
initial temperature of the system is obtained by solving the equation 



3.6 dN n 
ttR 2 a 4a cff (T l ) n dy 



Tf = ^TT^FTZ ~£ (13) 



where dN n /dy is the total pion multiplicity (~ 1.5 x charge multiplicity), Ra is the 
radius of the system, r, is the initial thermalization time and a e ff = (7r 2 /90) g c g(Ti). 
When the system is produced in the QGP phase, <7 c ff is replaced by gQGP which for 
two quark flavours is 37. If the quark and gluon masses are non-zero in the thermal 
bath then the effective degeneracy of the QGP, <7qgp, defined by Eq. ( [l2|) can be 
lower than 37, resulting in a higher value of Ti for a given multiplicity according to 
Eq. (|l3|). Note that the change in the expansion dynamics as well as the value of the 
initial temperature due to medium effects relevant for a hot hadronic gas also enters 
through the effective statistical degeneracy. The freeze-out temperature Tp is taken 
as 120 MeV. 



5 Results 

In Fig. |l|, the photon yield with QGP in the initial state is shown for three different 
values of the initial transverse velocity. The yield from hadronic matter during the 
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Figure 2: Total photon yield in Pb + Pb collisions at 158 A GeV at CERN SPS. The system 
is formed in the hadronic phase with the hadronic masses approaching zero at initial temperature 
T t = 200 MeV. 



mixed and hadronic phases are calculated with medium effects in the scaling scenario. 
All the three curves represent the sum of the thermal and the prompt photon con- 
tribution which includes possible finite fey effects of the parton distributions. The 
later, shown separately by the dot-dashed line also explains the scaled pp data from 
E704 experimented. We observe that the photon spectra for the initial velocity profile 
given by Eq. (|lTJ) with vq = 0.2 explains the WA98 data reasonably well. It is found 
that a substantial fraction of the photons come from mixed and hadronic phase. The 
contribution from the QGP phase is small because of the small life time of the QGP 
(~ 1 fm/c). 

To check whether the existence of the QGP phase essential to reproduce the WA98 
data, we have considered two possibilities: (a) pure hadronic model without medium- 
modifications, and (b) pure hadronic model with scaling hypothesis according to 
Eq.(|l]). In the former case, Tj is found to be ~ 250 MeV for r, = 1 fm/c and 
dN/dy — 700, which appears to be too high for the hadrons to survive. Therefore 
this possibility should be excluded. On the other hand, the second case with an 
assumption of Tj = T c (which is just for simplicity) leads to Tj ~ 200 MeV, at Tj = 1 
fm/c, which is not unrealistic. In this case, the hadronic system expands and cools 
and ultimately freezes out at Ty=120 MeV. The results for three values of the initial 
radial velocity including the prompt photon contribution are shown in Fig. |^. The 
experimental data are well reproduced for vanishing initial transverse velocity also. 
This indicates that a simple hadronic model is inadequate. Either substantial medium 
modifications of hadrons or the formation of QGP in the initial stages is necessary to 
reproduce the data. It is rather difficult to distinguish between the two at present. 

We will now compare the results of our calculation with the.-dilepton spectra 
obtained by the CERES/NA45 Collaboration. We have consideredlia dileptons with 
transverse momentum px above 200 MeV/c and an opening angle ee > 35 mrad. 
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Figure 3: Dilepton spectra for (N c h ) =250 for different initial states and mass variation scenarios. 



These are kinepatical cuts relevant for the CERES detector and are incorporated as 
described in E2l. In all the figures, the quantity (N c h) indicates the average number of 
charged particles per unit rapidity interval in the pseudorapidity range 2.1 < r\ < 2.65. 
In all the cases background due to hadron decays are added to the thermal yields. 

We display the results for dN c h/d/q — 250 in Fig. |j. In this case apart from the 
hadronic gas scenario we have also considered a deconfmcd initial state with Tj = 180 
MeV which evolves into a hadronic gas via a mixed phase. The observed enhancement 
of the dilepton yield around M ~ 0.3 — 0.6 GeV can be reproduced with the QGP 
initial state, once the variation of vector meson masses in the mixed and the hadronic 
phases are taken into account (solid line). The data is also reproduced by a hadronic 
initial state with Tj = 190 MeV in the universal mass variation scenario (dashed 
line). The p-peak in the dilepton spectra shifts towards lower M for universal scaling 
compared to QHD model, indicating a strong medium effect in the former case. 

The dilepton spectra for QGP and hadronic initial states for dN c h/d-q = 350 is 
shown in Fig. || Results with hadronic initial state and universal scaling of hadronic 
masses with temperature describes the data reasonably well. We see that with the 
temperature dependent mass from QHD model the low mass enhancement of the 
experimental yield cannot be reproduced. A good description of the data can be 
obtained by taking Tj = 200 MeV with QGP initial state for T c = 190 MeV. We also 
show the results due to large broadening of the p spectral function in the medium. 
The broadening of p is modelled by assuming the temperature dependent width as : 

r p (T) = r,(o)/(i-r 2 /T c 2 ). 

6 Summary and Discussions 

We have evaluated the high energy photon yield in Pb + Pb collisions at CERN 
SPS energies with two different initial conditions. In the first scenario, we start 
with the assumption of the formation of a QGP phase at an initial temperature of 
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Figure 4: Dilepton spectra for (A r c / l )=350 for different initial states and mass variation scenarios. 



~ 200 MeV and then the system continues through mixed phase and hadronic phase 
before freeze-out. In the second scenario, we assume a chirally symmetric phase 
where the hadronic masses approach zero at a temperature ~ 200 MeV and then the 
system evolves towards freeze-out. The effects of the variation of hadronic masses 
on the photon yield have been considered both in the cross section as well as in 
the equation of state. The photon spectra reported by the WA98 collaboration are 
well reproduced in both cases. We have also studied the dilepton yield measured by 
CERES experiment for two values of the charge multiplicity in Pb + Au interactions. 
It is observed that the data can be described by both QGP and hadronic initial 
states with an initial temperature ~ 200 MeV. We have assumed T£ = 1 fm/c at 
SPS energies, which may be considered as the lower limit of this quantity, because 
the transit time (the time taken by the nuclei to pass through each other in the CM 
system) is ~ 1 fm/c at SPS energies and the thermal system is assumed tabe formed 
after this time has elapsed. There are also uncertainties in the value of T c Eil, a value of 
T c ~ 200 MeV may be considered as an upper limit. Moreover, the photon emission 
rate from QGP given by Eq. (|), evaluated in Refs. BEl by resumming the hard 
thermal loops is strictly valid for g << 1 whereas the value of g is ~ 2 at T ~ 200 
MeV. At present it is not clear whether the rate in Eq. (0) is valid for such a large 
value of g. 
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